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ABSTRACT 

Aims. In recent years several in-depth investigations of the three Galactic Of?p stars were undertaken. These multiwavelength studies 
revealed the peculiar properties of these objects (in the X-rays as well as in the optical): magnetic fields, periodic line profile variations, 
recurrent photometric changes. However, many questions remain unsolved. 

Methods. To clarify some of the properties of the Of?p stars, we have continued their monitoring. A new XMM-Newton observation 
and two new optical datasets were obtained. 

Results. Additional information for the prototypical Of?p trio has been found. 

HD 108 has now reached its quiescent, minimum-emission state, for the first time in 50-60yrs. 

The echelle spectra of HD 148937 confirm the presence of the 7d variations in the B aimer lines and reveal similar periodic variations 
(though of lower amplitudes) in the He i A 5876 and He ii A 4686 lines, underlining its similarities with the other two prototypical Of?p 
stars. 

The new XMM-A'ewton observation of HD 191612 was taken at the same phase in the line modulation cycle but at a different orbital 
phase as previous data. It clearly shows that the X-ray emission of HD 191612 is modulated by the 538d period and not the orbital 
period of 1542d - it is thus not of colliding- wind origin and the phenomenon responsible for the optical changes appears also at work 
in the high-energy domain. There are however problems: our MHD simulations of the wind magnetic confinement predict both a 
harder X-ray flux of a much larger strength than what is observed (the modeled DEM peaks at 30-40 MK, whereas the observed one 
peaks at 2MK) and narrow lines (hot gas moving with velocities of 100-200km s"', whereras the observed FWHM is ~2000km s"'). 

Key words. X-rays: stars - Stars: early-type - Stars: individuals: HD 108,HD I48937,HD 191612 - Stars: emission-line 



Tjj- 1. Introduction IWalborn et al. L l2010h . Stars that appear otherwise normal but 

l/S^, . . . ., , displa'y strong C m A 4650 emission are now called "Ofc" stars 

The presence or C m /i 4o50A in enussion with a strength com- d'Walborn et al L 12010.) 

parable to the neighbouring Nm lilies is not a usual feature of O- jj^^^^ prototypical Of?p stars are HD 108, HD 148937, 

• star spectra. In the past, stars displaying this feature were clas- juT-iim^n t-u *u * » j- j u- * e ^ 

vo • .n , , i „, „ 1 nrf^ . f ^ ^ , „™ and HD 191612. They are the most studied ob ects of that cate- 

g , sified by | Walborn| (ll972l) m a separate category, dubbed Of?p. ^j^^^ ^^^^ questions still remain unanswered. 

In recent years, this category has attracted a lot of attention, mo j- i i » j u . . • r> i j n 

,, , ,■ ■ „ , -.^^o HD 108 displays correlated photometric, Balmer, and Hei 

, as exemplmed by numerous publications Walborn et al. , 2003, ,. ^-i • . • i. • u • • 

^-H . ^ . -' , , , ,, , • , line prohle variations: the star is brighter when the emission 

• • , 2004; Donatietal. , 2006; Howarthetal. , 2007; Hubri2etal., ,. ^ ^ , ^ . i p,nni - ^ \ Rnnou^ t-u 

onno M ' » 1 mm onn^ onm mno x -ru ^ a' "Uss are Stronger (Na ze et al. L 12001; NazeetalJ, g008b). The 
• 2008; .Naze et al. . ,2001, .2004. ,2007. ,2008a) . These studies un- 



X 



., , , ,. .J., ,j. variations appear recurrent, with a probable time scale of about 

veiled the peculiar properties of these stars (for a summary see , ^ , \ UnnA . t-u • • 

rrz — 7^ — rrrTr::^:T7n ii i ■ r , , ^ ■■ 55 years (Naze et al. U2OO60. The minimum emission state was 

INaze et all I2008bl) and led to a modmcation of the dehnition / ^ ... ^„„ . , ^ . i u„n/iK 

— T^Z'^r^nn u . . not vet reached in 2004 (Naz e et al. L 120041) . preventing a precise 

of these stars. The Of .'p phenomenon now covers stars present- / . ~. Tr 17 jj-.- i i. 

, . . , „ 1 TT ^ r.- evaluation of the recurrence time. Recently, an additional obser- 

ing recurrent spectral variations (in the Balmer, He I, C m. Si III , • u j • . r tit-; irvo j . c 

^ o vational result enriched the picture of HD 108: the detection of 

lines), strong C iiU4650A in emission (at least at some phases), ^ ^ magnetic field (observed strength is 100-150G, corre- 

narrow P Cygni or emission components in the Balmer hydro- jing to a potential dipolar field of 0.5-2 kG .Martins"^ 

gen lines and He i lines (at least at some ph ases), and U V wind |2Qlo|) 

lines weaker than for typical Of supergiants (jNazeetalJllQOSi behaviour of HD 191612 appears very similar to 

that of HD 108, but^with a shorter, well-defined period of 
537. 6d iHowarth et al. L l2007l and references therein). A strong 



* Based on observations collected at the Haute-Provence 
Observatory, at the La Silla and San Pedro Martir Observatories, 

and with XMM-Newton, an ESA Science Mission with instruments magnetic field (-220+38 G) was detected for this star by 

and contributions directly funded by ESA Member States and the USA Donati et al. | (|2006|), who further proposed an oblique magnetic 

(NASA). rotator model for the star. Following these authors, the long 

** Research Associate FRS-FNRS period would result from magnetic bra king of the stellar rota - 

*** Operated by the Association of Universities for Research in tion (for an analysis of such effects see lud-Doula et all |2009|) . 

Astronomy, Inc., under NASA contract NAS5-26555. Studies are now under way to check the exact geometry of the 
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magnetic field (G. Wade, private communication). Finally, it is 
interesting to note that this star has a lower-mass companion in 
a 1542d orbit (Howarth et al. , 2007). 

HD 148937 does not display the large photometric changes 
and line profile variations unlike the two other objects of this 
study. Instead, changes of limited amplitude (peak's amplitude 
varying by a few percents vs several 100% for the other proto- 
typical Of?p stars) were observed for the Ha line, with a pos- 
sible period of about 7d, during a long-term monitoring cam- 
paign dNaze et al. Ll2008a l). but the temporal sampling of the op- 
tical dataset was not optimized for finding such a short period 
and the resultant value could therefore be subject to large er- 
rors. In addition, no change was detected in the low-resolution 
data for other spectral lines, which differs from what is ob- 
served in HD 108 and HD 191612. Finally, a tenta tive detection 
of a m agnetic field was reported for this star by iHubrig et al.l 
(l2008l ; B, - -276 ± 88 G), but a spectro-polarimetric moni- 
toring throughout the 7d cycle is clearly needed to ascertain the 
strength and geometry of the field. 

In addition to dedicated optical campaigns, several XMM- 
Newton observations of these three objects were obtained. The 
derived properties are remarkably similar: overluminosities with 
respect to the 'can onical' Lx/Lbql- 10~^ relation, soft spectra, 
broad X-ray lines dNaze et al. "I I2004 120071 12008ah . Monitoring 
HD 191612 further revealed small variations of the X-ray flux. 
However, as all data were taken during the same orbital cycle and 
the same 538d cycle, the exact cause for these changes remained 
unclear 

Though the published studies have led to the characteriza- 
tion of many aspects of these Of?p stars, several questions still 
remain unanswered - usually, these questions are different for 
the different objects (e.g. on the variations of the X-rays of 
HD 191612, on the evolution of the optical spectra of HD 108, 
or on the variation timescale of HD 148937), due to the inho- 
mogeneous datasets reported up to now. In the past few years, 
we have continued our monitoring of these three peculiar ob- 
jects, with the objective of characterizing and better understand- 
ing their properties. This paper reports the results of this cam- 
paign. Section 2 presents the datasets and Section 3 the results; 
Section 4 summarizes the new information. 

2. Observations and data reduction 

2.1. Aurelie data 

Since 2004, when optical spectra of HD 108 were last published 
dNaze et al. LI2004I) . we have continued the monitoring of the star 
with the Aurelie spectrograph mounted on the 1 .52m telescope 
of the Haute-Provence Observatory (OHP, France). The detec- 
tor was a thin back-illuminated CCD composed of 2048x1024 
pixels of size 13.5 //m^. We used the grating # 3, reaching a re- 
solving power between 8000 and 10 000 in the blue and red part 
of the spectrum, respectively. With this instrumental setup, we 
obtained spectra in the 4450-4900A, 5480-5920A, and 6340- 
6780A ranges. Exposure times were typically of 20-40 min, re- 
sulting in a signal-to-noise ratio of about 200 on average. The 
data were reduced in a standard way using MIDAS. The wave- 
length calibration was done using ThAr exposures taken imme- 
diately before/after the stellar spectrum, and the spectra were 
normalized using low-degree polynomials. 

An additional echelle spectrum obtained in June 2009 at 
San Pedro Martir with the Espresso spectrograph was also made 
available to us by P. Eenens and L. Mahy. This dataset has been 
added to the analysis and discussion. 



2.2. Coralie echelle spectra 

Twenty echelle spectra of HD 148937 were obtained in 2008 
May 8-21 with the 1.2 m Euler Swiss telescope at La Silla 
(Chile), equipped with the Coralie spectrograph and an EEV 
2kx2k CCD with pixel size 15yum^. The Coralie instrument is 
an improved version of the Elodie spectrograph (iBaranne et al.L 
1996). These observations covered the spectral range 3850- 
6890A with a resolving power of 55 000. During the first week, 
two spectra per night were taken; during the second week, one 
spectrum per night was taken. The integration times were fixed 
to 30min (except for the first two exposures, which were each 
20 min long) and the typical S/N ratio at 5000A was 150. The 
data were first reduced with the Coralie pipeline and then nor- 
malized using low-degree polynomials. A telluric correction was 
applied for Ha and He 1/15876 using the telluric spectrum of 
iHinkle et all d2000l) and the task telluric within IRAF. 

2.3. XMM-Newton observation 

HD 191612 was observed anew by XMM-Newton on 2008 April 
2 (Rev. 1523, ObsID 050068201, PI Naze), for a total exposure 
time of 24 ks. The instrumental configuration, processing, and 
extraction regi ons were identical to the first four observations 
of the source dNaze et al. L 12007!) . We processed the data with 
the Science Analysis System (SAS) software, version 8.0; fur- 
ther analysis was performed using the FTOOLS tasks and the 
XSPEC software v 1 1.2.0. No high background episode (due to 
soft protons) was detected at high-energies during this exposure. 

In addition to EPIC data, XMM-Newton provides high- 
resolution spectra thanks to the RGS instruments. The RGS data 
were reduced with the SAS and then combined wit h the previ- 
ous observations following the procedure outlined in lNaze et al. I 
(2007). However, the improvement in signal-to-noise of the 
combined spectrum is not important as all observations are of 
similar duration. We therefore have no new results to report on 
the high-resolution spectrum of HD 191612, a quantitative leap 
in quality awaiting the advent of more sensitive observatories. 

3. New results on Of?p stars 

3.1. HD108 

The s pectrum of HD 108 has been followed intensively since 
1986 dNaze et al. L |200U |2004). While some (purely photo- 
spheric) lines such as Hen /1 4542 remain constant, the strength 
of Balmer, He: (e.g. at .1.14471,4713,6678), Sim (e.g. at /1 4552) 
and C III AA 4650 lines have decreased with time, with profiles 
changing from pure emissions or P Cygni profiles to strong ab- 
sorptions or greatly reduced emissions, as illustrated by Fig. [T] 
These lines thus appear as a blend of a constant photospheric 
absorption and a varying emission c omponent. A sim i lar c hange 
was observed about 55 years ago dNaze et al. L 1200 iL 1200 6) but 
due to the lack of detailed monitoring in the past - often, only a 
global idea of the profile (emission vs absorption) was recorded 
- the timescale of the variations could not be precisely assessed 
and some even questioned the recurrence of the phenomenon. It 
is therefore of utmost importance to see how the lines evolve and 
if quiescence is reached. 

Though the emission compon ent had nearly disappeared 
these last years (Naze et aO. 12004 V the emissions continued to 
decline until 2004. Now, HD 108 has finally reached its quies- 
cent state since the data taken between 2005 and 2 009 are very 
similar (Figs. \T\ and |2]i. Based on past behaviour dNaze et al. L 
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Fig. 1. Evolution of the line profiles in recent years for HD 108. 



1200 Ih . it is expected that the emission will strengthen again 
in the future, as is seen in HD 191612. We note however that, 
for HD 191612, the quiescent state occurs during one-third of 
the cycle, which would correspond to ~18yrs for HD 108. 
Assuming that the variations are due to an oblique magnetic 



rotator configuration, this quiescence timescale would however 
depend on the exact geometry of the system, which may be dif- 
ferent between the two objects despite their obvious similarities, 
from literature, it is only known that Hfi remained in absorption 
during about 10 years, while He 1/14471 did so during several 
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Table 1. Mean and dispersion of the measured RVs and EWs 
for HD 148937. 
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Fig. 2. Evolution as a function of time of a few quantities related 
to HD 108. Upper panel: EW (in A, in the range 4845^870A) of 
HjS. Middle panel: spectral type criterion log(|^^^). Lower 
panel: RVs (in km s" ' ) of He n 4542,4686. 



decades (iNaze et al. L 1200 ih . However, the exact evolution was 
not recorded and it cannot be concluded from the older data how 
long the lines remained constant, in a quiescent state. 

3.2. HD 148937 

Previously, a long-term monitoring (3 yrs, with monthly obser- 
vations) provided a large set of low-resolution spectra. These 
data revealed the va riability of 3 lines: Ha, Hj6, and Hen /1 4686 
dNaze et al. L l2008ah . However, this variability was of low am- 
plitude: the EWs measured for Ha, the most variable line, dis- 
play a dispersion which, though it was 10 times larger than that 
observed for the neighbouring narrow Diffuse Interstellar Band 
(DIB), correspond to a change of only 20% in EW (it also corre- 
sponds to a variation in the peak's amplitude of only 7% between 
the two extreme profiles), whereas the other prototypical Of?p 
stars have Ha lines varying from absorption to emission, with 
EWs and peak's amplitudes varying by >>100%. Moreover, the 
He I and C iii lines seemed relatively constant, a behaviour differ- 
ent from that of HD 108 and HD 191612 - though subde changes 
may have remained hidden in such low -resolution data. A period 
search further revealed a period of 7.03 1+0. 003d when study- 
ing closely the line profile variations of Ha. However, the data 
sampling, aimed at studying monthly-to-yearly variations, was 
not adequate for detecting such a short period, which thus re- 
quired confirmation with a more intense temporal sampling. A 
new, short-term monitoring was therefore undertaken, yielding 
20 high-resolution spectra over 2 weeks. 

For consistency, the radial velocities and equivalent widths 
(EWs) were estimated as in Naze et al. (2008a); they are given 
in Table [T] While the C iv and DIBs display constant EWs, the 
Balmer, He ii A 4686, and He i A 5876 lines show a larger disper- 
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0.416+0.038 


Ha 




-0.871+0.251 


DIB ^6613 


-3.3+0.5 


0.143+0.010 



sion. When plotted against time (Fig. O, obvious modulations 
are detected. Ha and He ii A 4686 lines display maximum emis- 
sion when Hj6, Hy, and He 1/15876 present a minimum absorp- 
tion. 

The variability detected when looking at EWs is co nfirmed 
by th e Temporal Variance spectrum (TVS, see Fullerton et al. L 
4996, and Fig.©, which compares, in each wavelength bin, the 
variations between spectra with the scatter expected from ran- 
dom noise. It is important to note that the variability has not the 
same amplitude for all lines of a given element: for the Balmer 
lines, the TVS is very significant for Ha, significant for Hfi, and 
barely significant for Hy. This might explain why the variability 
of, e.g.. He I A 447 1 is not detected - its amplitude, expected to be 
lower than for He i A 5876, is certainly hidden by the noise. Note 
that no significant variation of the C iii A 4650 lines is detected. 

Periodograms were ca lculated for t he 20- exposures dataset 
using the techniqu es of iHeck et al. I d 19851 see remarks in 
iGosset et al.1l2001l) in the wavelength interval where lines vary 
significantly (for details on the method see Rauw et al., 200li 
The results are shown in Fig.|5] the highest peak is located at a 
frequency of 0.1315d-i, 0.1353d-i, 0.1373d-', and 0.1397d-i 
for He lU 4686, Hfi, He U 5876, and Ha, respectively. These 
frequencies translate into a period of 7. 16-7.60±0.40d. The rel- 
atively large error, compared to the old dataset simply stems 
from the different observational timespans (2 weeks vs 2.5 yrs) 
while the period differences from one line to the other have two 
origins: the noise, which affects more the fainter lines such as 
He I /1 5876, and the different amplitudes of the variations (the 
signal is indeed more difficult to catch in case of low ampli- 
tude changes, e.g. for He 1/15876). The Ha line, which is both 
strong and the most variable line, yields the most reliable pe- 
riod, i.e. 7. 16±0.40d from the sole new dataset. Taki ng the errors 
into acc ounfl this agrees with the results reported bv Naz e et al. I 
(l2008ah . The 7d period was thus not an artefact from our inad- 
equate time sampling. HD 148937 therefore appears quite sim- 
ilar to the two other Galactic Of?p stars, though with a shorter 
period, a smaller amplitude of the variations and with the sole 
exception of the apparent constancy in the C iii lines. A spec- 
tropolarimetric monitoring should be undertaken to see if this 
could be related to e.g. a low inchnation of the star's rotation 
axis. 



' Note that combining the two datasets is not possible due to the 
very different spectral resolutions (i?=2300 vs 55000). As the best 
period estimate co mes from the longes t time coverage, the value of 
7.031+0.003d from iNaze et al. I ( l2008ij) should be considered as the 
best measurement on an absolute scale. 



Yael Naze et al.: New findings on the prototypical Of?p stars 



5 



Strong emission 



• Ha 

o HelI4686 



cP o 



small emission 



4595 4600 4605 

HJD - 2 450 000 



strong absorption 



H7 



small absorption 



HIS 



HeI5876 



CIV5812 



DTB 6613 

I I I I ^ 

4600 

HJD - 2 450 000 



o 8 



Fig. 3. Evolution with time of the EWs (expressed in A) for some lines of HD 148937. Left: The strongest emission lines (Balmer 
line of Ha and He ii A 4686) ; the emission strength increases from bottom to top. Right; Some absorption lines ; the absorption 
increases from bottom to top. The He i A 5876 and Ha data were corrected for the contamination by telluric lines. 
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Fig. 4. Left: Average spectra and TVS of He ii A 4686, H^S, He U 5876, and Ha for HD 148937. The He U 5876 and Ha data were 
corrected for the contamination by telluric lines. The horizontal dotted line in the lower panels represents the 99% confidence level. 
Right: The 20 Coralie observations of the Ha line (from May 8 to May 21, note the high noise in the data of May 19, which reflects 
in the EW measurements) and the extreme states of Ha, He i A 5876, and H/3. The black solid lines correspond to the spectra of May 
09 (max emission state) and May 12 (min), the dotted red lines to those of May 15 (max) and May 20 (min). The colored version of 
the figure is available on-line. 



3.3. HD 191612 
3.3.1. X-ray data 

The previous observations of HD 191612 (iNaze et aini2007l) re- 
vealed changes in the X-ray flux but the origin of this variability 
was unclear. To get a better understanding of this phenomenon. 



we undertook a new XMM-Newton observation, whose schedul- 
ing was carefully chosen. While the previous observations were 
taken at orbital phases 0ori=O.83-O.96 and at phases 0c,.c=O.O9- 
0.44 in the 538d cycle, the new exposure corresponds to a phase 
0cvc=O.13 in the line-profile/photometric cycle of 538d and a 
phase (^or<)=0.55 in the orbital cycle (Fig. |6l ephemeris taken 
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Fig. 5. Left, from top to bottom: Raw periodogram of Ha as observed during the May 2008 campaign, periodogram after prewhiten- 
ing by the period of 7. 16d, and associated spectral window in the case of HD 148937. The latter is directly related to the sampling 
of the time series. Right, from top to bottom: Periodograms of He ii A 4686, H/?, and He i A 5876. The power spectra were calculated 
and averaged in the interval of significant variability (see Fig. |4] and top label of each panel). The He 1/15876 and Ha data were 
corrected for the contamination by telluric lines. 



from lHowarth et al. l2007h . The new data thus enable us to check 
whether the flux variations are linked to the circumstellar, pos- 
sibly magnetic, structures surrounding the star or to colliding 
winds in a binary. In the former case, the new data should be 
nearly identical to the previous observation taken at a similar 
phase (pcyc- On the contrary, in the latter case, a dramatic varia- 
tion is expected. In a colliding wind binary, variations of the flux, 
correlated with orbital phase, are expected to arise from changes 
in the separation (in eccentric systems) and/or in the absorb- 
ing column (when the winds are different). Indeed, such vari- 
ations have been observed in several systems (for a review, see 
[Oudel & Naze, 2009). In HD 191612, the previous observations, 
taken before periastron, revealed a 20% change in reddening- 
corrected flux (see Table |2]l for a phase difference Acporh of only 
0.12, corresponding to a change in separation of 36%. The new 
observation was taken two-thirds of the orbit later, near apastron: 
the change in separation compared to previous data has thus in- 
creased (up to 60%), and so should the difference in flux, if truly 
linked to the colliding wind emission. 

The new data were processed following the same procedure 
as for the previous data, but using more recent XMM-Newton re- 
duction software. The count rates derived from a run of the task 
edetectjchain amount to 0.122±0.002cts s"' for EPIC MOSl, 
0.119+0.002ctss-' for EPIC MOS2, and 0.392±0.005cts s"' 
for EPIC pn. The evolution of the count rates and hardness ratio 
is shown in Fig. |6l The new data agree very well with previous 
observations taken at the same phase of the 538 d cycle, both in 
strength (similar count rates) and in hardness (similar hardness 
ratios) ; they are in addition clearly different from the results 
of previous observations obtained at the latest phase (where the 
count rates are ~35% smaller, corresponding to a ~20cr differ- 



ence, and the hardness ratio HRi measured on pn data is lower 
by 7cr). 

The best spectral fit is further provided in Table |2] 
together wi th the previous results for comparison (repro- 
duced from iNaze et al. 1 l2007h . The fitted model has the form 
wabs (A^fjj) * [wabs (A^f ) *mekal ( kri)+wabs(A^») '-mekal (k^a) ] , 
with wabs(A^j^,)= 3.4 x 10^' cm^^ jp ipias & Savage , 1994) and 
using solar abundances for both absorptions and o ptically-thin 
plasma emissions - i.e. the same model as used in iNaze et al. I 
(2007) for homogeneity reasonfl Quoted fluxes are in the 
0.4-1 0.0 keV energy range, and the unabsorbed fluxes /™*'* are 
corrected only for the interstellar absorbing column. For each 
parameter, the lower and upper limits of the 90% confidence 
interval (derived from the error command under XSPEC) are 
noted as indices and exponents, respectively. The normalisation 
factors are defined as J rieriudV, where D, rig and riu 
are respectively the distance to the source, the electron and 
proton density of the emitting plasma. The phases <porb and 
(pcyc correspond to the phases in the orbital cycle and the 
53 8 d line profile/photo metric cycle, respectively (ephemeris 
of iHowarfli et al. I l2007l) . Again, the spectral fits to the 2008 
data agree very well with the previous observations taken at 
the same phase of the 538 d cycle, yielding the same X-ray 
luminosity, 9x10^^ ergs"' (for a distance of 2.29 kpc), and same 
\og[Lx/LBOL], -6.1. 

The conclusions are obvious: despite the very different or- 
bital phases, the new data are nearly identical to those taken 3 
years before at a similar phase in the 538d cycle. The modula- 

- As in lNaze et al. I l l2007h . we noted the two sets of parameters giving 
a good fit, but we recall that a fit by a differential emission measure 
(DEM) model clearly favors the 'cool' fit. 
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tion of the X-ray emission thus appears directly correlated with 
the optical line-profile variations, a situation very similar to that 
observed for 6^ Ori C (,Gagne et al. .2005; .Naze et al. ,2008b.) . 

3.3.2. MHD modelling 

To better understand the properties of HD 191612, a fully- 
dynamic numerical magneto-hydrodynamic (MHD) model, tak- 
ing into account the impact on the wind of the magnetic field 
detection, is needed. We used the publically available ZEUS- 
3D code (Stone & Norman, 1992) to evolve a consistent dynam- 
ical solution for a line-driven stellar wind from the star with 
a dipolar surface fi eld. The stellar par ameters of HD 191612, 
adopted mostly from lNaze et al ] (l2008bl and references therein), 
are listed in Table |3] 

HD 191612 is a rather slow rotator w i th a surface rotatio n 
speed <50km s ' (iHowarth et al. L 120071 iNaze et all l2008bh . 
Such a slow rotation is unlikely to have any significant dynami- 
cal effects on the stellar wind. Therefore, the wind was assumed 
to be azimuthally symmetric with the magnetic field axis aligned 
with the rotation axi^. This is essentially a '2.5-D' formulation 
allowing for non-zero azimuthal components of both the mag- 
netic field and velocity v^, while still assuming that all quan- 
tities are constant in the azimuthal coordinate angle <p. 

Much of the numerica l procedures in this simulation fol- 
lows |ud^Douiaeta (120081) . The models thus include rotational 
effects but, for the first time, we also use the full energy 
equation, with adiabatic index y=5/3 and radiative cooling of 
[MacDonald & Bailev ( 1981, which assumed solar abundancefl 
At the initial time f = 0, the wind is represented by a relaxed 
non-magnetic and non-rotating model, but we then simultane- 
ously introduce a dipole magnetic field and a surface rotation 
at the lower boundary, both defined relative to a common polar 
axis. The wind and magnetic field are then let free to compete 
with each other. To ensure that there are no numerical artefacts 
arising from the initial conditions, we run our model for about 
1.5 Msec, which is at least 50 fluid crossing times. 

Immediately after f = 0, the magnetic field channels wind 
material towards the tops of closed loops near the equator. There, 
the collision with the opposite stream leads to strong shocks, like 
in "magnetically confined wind shocks" (MCWS) presented by 
[Babel & Montmerle (1997). Denser material then cools radia- 
tively very quickly leading to a dense disk-like structure (Fig. 
|7]i. Some of the less dense post-shock material within the closed 
magetic loops remain very hot (ca. 80MK) for a long period of 
time. In this model, the Alfven radius or magnetosphere extends 
about 3-4 stellar radii above the stellar surface. 

Initially mass builds up in the equatorial region below the 
Keplerian radius estimated to be at about 4.5/?,, but then there 
appear repeated episodes of infall of inner disk material back 
onto the stellar surface. Over a somewhat longer timescale, there 
appears another, somewhat different kind of disruption, one that 
starts higher up, closer to the Alfven radius. This is characterized 
by outward ejections of the upper disk mass, leading to 'centrifu- 
gally driv en m ass ejections'. Such ejections can lead to X-ray 
flares (see lud-Doula et al.ll2006l) but these short-term variations 
cannot be detected in the observations of HD 191612, since they 



' This simplifies the computation and makes no assumption on the 
true geometry of the system, i.e. an oblique rotator model is not ex- 
cluded ; it can easily be simulated here by changing the viewing angle 
on the system. 

Note that models used bv lGagne et alj ( 120051) made use of the full 
energy equation but did not include rotation. 



Table 3. Stellar parameters used for the MHD simu lation of 
HD 191612, adopted mostly from iNaze et"an (l2008bl and ref- 
erences therein). 



Parameter 


Value 


R 


1.015x10 cm = 14.4 Rq 


M 


7.0xlO^^*g = 35Mo 


L 


7.697xl03*ergs-' =2x105Lo 


a 


0.55 


Q 


2000 


6 


0.0 


M 


1x10-* Mo yr-' 




50km S-' 




1500 G 




31 ~ VlOOO 



lack the required sensitivity to analyze the lightcurve on very 
short timescales. 

Fig.|8]shows the differential emission measure (DEM, or the 
relative amount of material in a given temperature bin) calcu- 
lated in the MHD simulation. The DEM peaks at very high tem- 
peratures (about lOMK and larger), suggesting that the very hot 
gas dominates the X-ray emission of the system. Of course, the 
material in the unshocked regions remain at the effective temper- 
ature of the star due to the intense UV radiation from the parent 
star that keeps the wind nearly isothermal. 

3.3.3. Comparison of models and observations 

Applying a magnetically-confined wind model to HD 191612 
yields four clues. First, the X-ray emission should clearly be 
dominated by very hot plasma, as was the case for 0' Ori C (see 
Fig. 6 of iGagne et al.L l2005h . Second, on average, the total 
EM calculated for HD191612 in the 1-lOOMK range amounts 
to 5.3x10^^ cm-^. This agrees well with the MHD simulation 
of 0' Ori C, considering the difference in the stellar parameters 
and modelling. Indeed, the simulated EM for this star, in t he 
same temperature range, reached 8.2x10^^ cm"^ dOagne et all 
2005), i.e. a factor 6.5 lower than for HD 191612 which is close 
to the factor of 4 expected from the EMs scaling with the third 
power of the radius. Considering an average cooling factor MJ) 
of ~ 2 X lO"^-' erg cm-' s"' (Raymond et al., 1976), this yields a 
totally unabsorbed X-iay luminosity of ~ 8 x lO-''' erg s"' in the 
0.5-10. keV range. Third, in the MHD models, the hot gas gen- 
erally moves with velocities of about 100-200km s"' and only 
very few hot plasma reach 400km s ', whatever the viewing 
angle on the equatorially-confined region. This is unsurprising 
since the hot gas always remains close to the star: the associated 
X-ray lines are thus expected to be narrow, again a similar result 
as that found for 6*' Ori C. Finally, the X-ray emission should be 
modulated throughout the cycle of 538d, as the viewing angle on 
the disk-like structure changes. 

The observations (both old and new) do not fully agree with 
this picture. On the one hand, the emission is modulated with 
the 538d cycle, as expected for MCWS. The observed X-ray 
emission is indeed bright (Lx-1-9xlQ^^ erg s"' when corrected 
for the interstellar absorption), in fact on the luminous edge 
compared to "normal" O-type stars (for O stars observed with 
XMM-Newton, log [Lx/^-bol] =-6 .45 with a dispersion of 0.51 
dex, see Naze 2009). Considering the uncertainty in the spec- 
tral models (see Table |2]i, the observations yield a totally un- 
absorbedX-ray luminosity of ~8-80xl0^^ergs-', for the 'hot' 
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Table 2. Best-fitting models and X-ray fluxes at Earth for each XMM-Newton observation of HD 191612. 



Date 


Rev. 




(porb 






norm I 




kr. 


norm2 




Jx 


xunabs 
.'X 










lO'''^ cm ^ 


keV 


10 ''cm ^ 


cm ^ 


keV 


I (-V 3 5 

10 ""cm ^ 




(W-^ erj 


;cm s ) 


'cool' model 
























05/04/05 


975 


0.09 


0.84 


0.5 IJJ;^^ 




7.07i5,^ 




1.27!j« 


1 nci 1 7 1 


1.16 (411) 


7.0 


14.1 


17/04/05 


981 


0.12 


0.84 


^7^'^" 
0.25 


27 


--1 1 c-"^ SI 

2 1 S-^-°^ 


V 1 


1 r\ t 1 07 

1 01 




1.10 (681) 


6.8 


14.3 


02/06/05 


1004 


0.20 


0.87 


() 4i 


25 


^•'°i:78 




1 22 


0.89i:^i 


1.15 (460) 


6.0 


12.5 


08/10/05 


1068 


0.44 


0.96 


61 


20 


11 ^18 1 

1 3 6 


1 1 fi8 

'-■■'^1.14 


1 -601:15 


() 66 

"•""o.59 


1.31 (555) 


5.3 


11.0 


03/04/08 


1523 


U.l J 


U.JJ 


n 1 56 

0.51o4g 


n 740-25 

0.24^23 


6.67^,5 


1 1.41 




1.06(,,7 


l.jL {jyL ) 


n 7 


1 A A 
14.4 


'hot' model 
























05/04/05 


975 


0.09 


0.84 


qO.05 
"•0. 


64" '''' 

"■'-'^0.62 


0.28J5;27 


Q 0.03 
"•0. 


2.35|« 


51"-^^ 

"•-"^0.48 


1.19 (411) 


7.6 


14.4 


17/04/05 


981 


0.12 


0.84 


"•0. 


60" ''' 




A 1 0.87 




95''5 

"•=^-^0.86 


1.08 (681) 


7.0 


13.9 


02/06/05 


1004 


0.20 


0.87 


"■0. 


51"" 

"■-^^0.47 




0.82^;?i 


1 07 '-2' 
^■"'i.oi 


92'"' 


1.16(460) 


5.8 


12.0 


08/10/05 


1068 


0.44 


0.96 


qO.OI 
"•0. 


61"" 

"•"^0.60 


23"-2'* 

2'' 


0-r' 


2.532 3g 




1.39 (555) 


5.4 


11.0 


03/04/08 


1523 


0.13 


0.55 




61" ^^ 

"•"^0.60 


33"-^^ 

0.32 


o.57g:S 


1 64'-'* 

^•"^1.53 


73"-^'' 

"• '-'0.66 


1.31 (591) 


7.3 


14.1 




Fig. 6. Left: Position of the new XM M-A^ewfow observation in the cycles of HD 191612. The orbital cycle, calculated using orbital 
parameters of lHowarth et al. I ( l2007h and with the primary (resp. secondary) curve drawn in solid (resp. dotte d) line, is shown on top 
while the bottom panel shows the line profile variation cycle of 538d, from the analytical approximation of lHowarth et al. I (l2007h . 
The phases of the previous observations are shown by the thin short-dashed lines and the phase of the new observation by the thick 
long-dashed line. Right: Evolution of the count rate in the 0.4-10. OkeV band (top) and of two hardness ratios (middle, bottom) with 
phase from the 538d cycle. Open symbols refer to the old observations, filled red symbols to the new dataset; triangles and squares 
were used for the EPIC MOS data (note that the MOS count rates are multiplied by a factor of 3 in this figure) and circles for the 
EPICpn data. The hardness ratio HRx and///?2 are defined as {M-S)I{M + S) and {H - M)I{H + M), respectively, where S is the 
count rate in the 0.4-1. OkeV band, M in l-2keV, and H in 2-lOkeV. The colored version of the figure is available on-line. 



and 'cool' model respectively. Total emission meas ures, derived 
from t he normalization factors of Table |2] or as in iGagne et alj 
(l2005h for a spectral fitting considering a single absorbing col- 
umn for both thermal components, amount to 3-7 x 1 0^* cm"^ for 
an equivalent 'cool' model and ten times lower for an equivalent 
'hot' model. Theoretical predictions thus agree rather well with 
the 'cool' model, which is also the one favored by the DEM 
modelling dNaze et al. L l2007l) . 



On the other hand, however, the X-ray emission of 
HD 191612 displa ys broad fines (FWHMohs ~2000km s ^, 
iNaze et al H2007h " and is far from being as hard as predicted 
by the model. In fact, a fit of the EPIC data by a DEM model 
yields a strong emission peak at about 2 M K, not beyond 10 MK, 
as shown for example in Fig. 3 of ,Naze et al. I (l2007h . In ad- 
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Fig. 7. Snapshots of density (left) and temperature (right), in cgs units on a logarithmic color-scale, shown with field lines (solid 
lines) at arbitrarily chosen time t=1300ks for the 2D MHD simulation of HD 191612. the x and y axes are labelled in units of 
the stellar radius. Note the unsteady nature of the dense equatorial region which undergoes episodic emptyings and refills. Within 
the magneto sphere (i.e. the closed loop region), the opacity at lOA or 1.2 keV may reach high values (t=10 for /(■=30cm^ g"', see 
ICohenll2009l a density of 10 g cm and a disk height of 0.03 R,) but it is totally negligible in the dense equatorial region outside 
magnetosphere. The observations show only a small absorption in addition to the interstellar extinction and it is not significantly 
varying with phase (i.e. there is no clear evidence of excess absorption in t he magnetic equato rial plane, or the effect is much smaller 
than would be expected from a dense equatorial cooling disk, as found in lGagne et alj|2005l) . 



Fig. 8. Volume emission-measure distribution (DEM) per log T= 0.1 bin for our MHD simulations of HD 191612. The MHD 
simulation is used to compute the emission measure per unit volume per logarithmic temperature bin, then integrated over three- 
dimensional space assuming azimuthal symmetry about the magnetic dipole axis. Left: DEM as a function of time and temperature 
; the periodic variations due to build-ups and infalls are clearly seen. The colored version of this figure is available on-line. Right; 
Time-averaged DEM as a function of temperature. 



dition, the average temperatur^ derived in the 2XMM survey 
for HD 191612 amounts to ~0.3 keV and is thus quite quite nor- 
mal (85% of the stars in the 2XMM have < kT > below 1 keV, 
[Naze 2009) - it does not resemble the extreme case of 6*' Ori C 
for which < kT > amounts to 2.5 keV in the same 2XMM sur- 
vey and for which the MCWS model works well. The relatively 
soft character of that emission despite the large magnetic con- 
finement parameter might be reminiscent of the characteristics 
observed in some B-type stars (Giidel&Naze, 2009): further 
investigation of the MCWS model and magnetic field geome- 
try is thus required to understand the high-energy properties of 
the Of?p stars. 



4. Conclusions 

Our continued monitoring in the X-ray and optical domains of 
the prototypical Of?p stars has now clarified some of their ob- 
servational properties. 

For HD 108, the Balmer, C iii .14650, He i, and Si iii lines dis- 
play minimum emissions since 2005: the star has thus finally 
reached its quiescent state, for the first time in 50-60yrs. Long- 
term monitoring of this quiescence will further constrain the ge- 
ometry of the line emission region. 

For HD 148937, the presence of the 7d variations in the 
Balmer lines is confirmed by a short-term monitoring. During 
this cycle, the amplitude of the emission peak in Ha, the most 
variable line, changes by only 7%. Thanks to the high spectral 
resolution of the new data, similar periodic variations, though 
of even smaller amplitude, are detected in the He i /1 5876 and 
He 11/14686 lines. These detections underline the similarities of 
the star with the other two prototypical Of?p objects. 

For HD 191612, a new XMM-Newton observation clearly 
shows that its X-ray emission is modulated by the 538d period 
and not the orbital period - it is thus not of colliding-wind ori- 
gin: the phenomenon responsible for the optical changes appears 
also at work in the high-energy domain. Contrary to observa- 
tions, our MHD simulations predict a hard spectrum, dominated 
by very hot plasma (>10MK) and narrow lines. Further refine- 
ments to the modelling will be needed in order to fully explain 
this discrepancy. 



' 2 kTj X nornij/ 2 nornij for models using a single absorption in 
front of the thermal emission components 
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